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The finite size and surface roughness effects on the magnetization of NiO nanoparticles 
is investigated. A large magnetic moment arises for an antiferromagnetic nanoparticlc 
due to these effects. The magnetic moment without the surface roughness has a non- 
monotonic and oscillatory dependence on R, the size of the particles, with the amplitude 
of the fluctuations varying linearly with R. The geometry of the particle also matters 
a lot in the calculation of the net magnetic moment. An oblate spheroid shape particle 
shows an increase in net magnetic moment by increasing oblateness of the particle. 
However, the magnetic moment values thus calculated are very small compared to the 
experimental values for various sizes, indicating that the bulk antiferromagnetic structure 
may not hold near the surface. We incorporate the surface roughness in two different 
ways; an ordered surface with surface spins inside a surface roughness shell aligned due 
to an internal field, and a disordered surface with randomly oriented spins inside surface 
roughness shell. Taking a variational approach we find that the core interaction strength 
is modified for nontrivial values of A which is a signature of multi-sublattice ordering for 
nanoparticles. The surface roughness scale A is also showing size dependent fluctuations, 
with an envelope decay A ~ i?" 1 / 5 . The net magnetic moment values calculated using 
spheroidal shape and ordered surface are close to the experimental values for different 
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1. Introduction 

Antiferromagnetic nanoparticles have been receiving a refreshed research attention 
over the last few years. These are considered as better candidates for exhibiting 
the magnetization reversal by quantum tunneling 1, due to their small magnetic 
moment as compared to the ferromagnetic nanoparticles. The magnetic proper- 
ties of nanoparticles are dominated by fini te-siz e effects, and the surface anomalies 
such as surface anisotropy and roughness EGOS As the particle size decreases, the 
fraction of the spins lying on the surface of a nanoparticle increases, thus, mak- 
ing the surface play an important role. The reduced coordination of the surface 
spins causes a symmetry lowering locally, and leads to a surface anisotropy, that 
starts dominating as the particle size decreases. Thus, an enhancement of surface 
and interface effects make the antiferromagnetic nanoparticles an interesting area 
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Nickel Oxide (NiO) has been considered as a prototype for antiferromagnetism, 
as it is one of the first few materials in which antiferromagnetism was 

studiedPlOne 

of the first serious concerns with NiO nan opa rticle was, evidenced from the experi- 
mental study of Richardson and MilliganpSl that these nanoparticles show a large 
magnetic moment as the size becomes smaller than lOOnm, apart from anomalous 
behavior of the magnetic susceptibility. It was also found that the exchange coupling 
between the surface spins and the antiferromagnetic core spins causes an exchange 
bias phenomenon in these finite-sized particles. This phenomenon is responsible for 
the observed shifted hysteresis loop after field cooling in NiO nanoparticles P This 
interface effect is very much size dependent. A large loop shift ( > lOKOe) and co- 

ercivities at low tempe rature has been reported for the intermediate sized particles 
(22nm-31nm).12ini 

Winkler et al ^ reported that for 3nm particles the magnetization curves are 
reversible above T ~ 170K, but a hysteresis behavior is observed at lower temper- 
atures. According to their observation, a large surface anisotropy is responsible for 
the anomalies in the shape of the hysteresis loop at low temperatures. They found 
that with a decrease in the temperature, a progressive blocking of the core particle 
moments starts off, and it is followed by a growth of spin clusters at the particle's 
surface below 40K, and finally their collective freezing in a cluster glass-like state 
at 15K. 

The net magnetic moment of antiferromagnetic nanoparticles has been a subject 
of research interest from a long time. Neel in 1961 suggested that fine particles of 
antiferromagnetic materials exhibit weak ferromagnetism and superparamagnetism. 
He argued that the permanent magnetic moment in these antiferromagnetic fine 
particles is due to incomplete magnetic compensation between the atoms on the 
two sublattices 'A' and 'B', which are identical in every respect, except that the 
atomic moments in B sublattice are antiparallel to that in A sublattice. Neel con- 
sidered three general cases as shown in Ref. If the uncompensation of spins 
occurs randomly in a particle, then the number of uncompensated spins p will vary 
as p ~ Ti.2, where n is the number of spins. If the spins are arranged in such a 
way that the ordered structure consists of odd number of ferromagnetic planes of 
A and B atoms, then p ~ ns . Finally, if each plane consists of equal numbers of 
A and B atoms and the structure consists of incomplete top and bottom planes, 
then we would have p ~ ni . Richardson et al^^ showed that p ~ ns from the size 
dependence of susceptibility in NiO nanoparticles. Thus, according to the Neel's 
model, the magnetic moment /x for NiO nanoparticles var ies as [i ~ n3/ijvi 2 +- Weak 
ferromagnetism was later confirmed by experiments on fine particles of NiO. 
For extremely fine particles they reported the behavior to be superparamagnetic. 
However, Tiwari et al® argued that the NiO nanoparticles behave like a superspin 
glass, which is attributed to a surface spin disorder. Some authors accredited the 
large magnetic moment in NiO nanoparticle to nonstoichiometry, an existence of 
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small superparamagnetic metallic nickel clusters in NiO part icle or the presence of 
Ni 3+ ions within the NiO lattice PS However Richardson et aft^ confirmed that the 
presence of Ni 3+ in NiO do not contribute significantly to the magnetic moment of 
NiO nanoparticles. More recently Yi et al ^ investigated the size dependent mag- 
netic properties of NiO nanostructures using experimental and first principle study. 
They reported that NiO clusters with a size upto lnm indicate ferromagneticlike 
interactions with high magnetizations, and NiO nanocrystals with a particle size 
over 2nm possess uncompensated magnetization. 

The magnetic moment per particle for NiO has been investigated experimentally 
by Kodama et From extrapolation of hK magnetization curves from a large 
magnetic field to zero field, they found 700 ub per particle for particles of size 
15nm, while the Neel's two-sublattice model — predicts a magnetic moment of 
about 80 fiB- For the particles of size 3nm, Winkler et al ® experimentally found 
the magnetic moment per particle to be 500 /is, whereas for this particle size, the 
Neel's model predicts a magnetic moment of 20 \j,b- 

This discrepancy, between the magnetic moments experimentally observed and 
those predicted by the Neel's two-sublattice model has been a serious question from 
a long time. Kodama et al® have shown from numerical modeling that a reduced 
symmetry on the surface of the nanoparticle actually causes a fundamental change 
in the magnetic order which results in a multi-sublattice structure. Monte Carlo 
studies for antiferromagnetic nanoparticles by Zianni et also reveals a distinct 
magnetic role of surface and core spins. Recently it has been pointed out that the 
roughness at the surface layer gives r ise to higher magnetic response for the surface 
spins than the core spinsP2E2223E5] 

In view of these studies, we investigate the large magnetic moment in NiO 
nanoparticle by invoking a different ordering for surface spins than the bulk Neel 
state ordering for core spins. 

The outline of the present manuscript is as follows: In Sect.[2]we discuss a model 
for bulk NiO, followed by a discussion on the finite-size effect in the magnetization 
of NiO nanoparticles using spherical as well as spheroidal geometries in Sect. [3] We 
will discuss surface effects and ordering of surface spins beyond Neel state ordering 
for nanoparticles in Sect. [4] and Sect. [5] is devoted to conclusions. 



2. The model 

The crystal structure of bulk NiO has been comprehensively investigated in the 
literature using x-ray diffraction methodESEll. It has been found to be face centered 
cubic (fee), with a Neel temperature of 523K. Each Ni atom has twelve nearest 
neighbors and six next-nearest neighbors. The lattice parameter has been found ^Sl 
to be 4.1758A at 297K and 4.17051 atT^OK 

The magnetic structure of NiO has been well establish ed to be fcc-II by the 
work of Shull et al — and further by Roth et a /] 29 | 30 | 3 1] a ^ orn | c S p ms are 
stacked ferromagnetically in (111) plane but aligned antiferromagnetically in (111) 
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Fig. 1. The magnetic configuration of bulk NiO where closed-packed ferromagnetic sheets of 
spins are stacked antiferromagnetically along the direction perpendicular to the sheet which is 
(111) direction for bulk NiO. The yellow (white) spheres denote up spins and the blue (black) 
spheres denote down spins. 



directions. The direction of alignment of the spin moments has been found to be 
(112) directions. The magnetic configuration of bulk NiO is shown in Fig. [I] 

The neutron diffraction studies by Hutchings et al E2J confirmed that the predom- 
inant interaction in NiO is a large next-nearest neighbor antiferromagnetic exchange 
interaction J nnn = 221K linked by 180° superexchange path Ni 2+ — 2 ~ — Ni 2+ . 
The nearest-neighbor interaction is linked by 90° path Ni 2+ — 2 ~ — Ni 2+ , which is 
much smaller in strength. Due to the lattice contraction, there is a slight difference in 
the exchange interaction between the nearest neighbors in the plane, J~ n = 16.1-KT, 
and between nearest neighbors out of the plane, j£ n = \b.7K. Hutchings et al 
also used an orthorhombic form for the anisotropy Ef = Ki(si.x) 2 + Ki(si.z) 2 with 
x is the easy axis direction (112) and z is the hard axis direction (111), where Si 
is the atomic spin at site i. The anisotropy constants are gives as K\ = 1.13K and 
Ki = .06K. Since Ki is much smaller than K±, we use K\ as anisotroy constant 
and z as anisotropy axis. 

The spins in the NiO interact via Heisenberg exchange interaction. The Hamil- 
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Fig. 2. A portion of (111) plane of the bulk NiO is shown. A circle, whose centre is one of the 
lattice points itself is drawn. The circle shown is a cross section of NiO nanoparticle. The lattice 
points on the perimeter of the circle are highlighted. Some of the lattice points are lying just inside 
the perimeter and some just outside. These lattice points are responsible for the fluctuations. 

tonian of the system in the presence of an external magnetic filed H is 

7~L = Jnnn Si.Sj ~ J nn Si.Sj + J^ n Si.Sj 

(ij) (ij) <W) 

-K^isaf-H-Y.^- (!) 

i i 

The first term represents the dominant antiferromagnetic next-nearest neighbor ex- 
change energy. These next-nearest neighbors lie in the adjacent planes just above 
and just below the plane consisting the spin Sj, e.g. for each spin in the plane B in 
Fig. [TJ three of the six next-nearest neighbors lie in plane A while other three lie in 
plane C. The second term represents the ferromagnetic nearest neighbor exchange 
energy which determines the interaction of Si with six of the twelve nearest neigh- 
boring spins lying in the same plane as spin Sj. The third term is antiferromagnetic 
nearest neighbor interaction energy which represents the antiferromagnetic inter- 
action of Si with six nearest neighbors lying in the planes other than the plane 
containing spin Sj, e.g. each spin in the plane B, has three nearest neighbors in 
plane A and three in plane C. The fourth term represents the uniaxial anisotropy 
energy and the last term is the Zeeman energy. 

The most dominant term in the above Hamiltonian is the first term which supports 
antiferromagnetic order. Thus, in the bulk we have a Neel state ordering where 
spins are stacked ferromagnetically in (111) plane but aligned antiferromagnetically 
in (111) directions. Though the Neel state in the bulk has zero magnetization, how- 
ever, as the size of the particle becomes smaller, the Neel state ordering shows a 
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Fig. 3. The fluctuation in the total number of spins in a sphere of radius R is shown. The solid 
line shows the fluctuation A^^ere obtained from analytical expression Eq. JHJ for the terms up to 
G = 400 and Eq. J7| , while the dotted line is the exact numerical counting of the fluctuation from 
Eq. 10 and Eq. (Q. 

magnetization, since the spins lying near the surface do not cancel out. We will 
examine below the finite size effect on the Neel state magnetization. 

3. Finite-size effects in the Neel state 

We consider various geometries for antifcrromagnctic NiO nanoparticles. The crystal 
structure of NiO nanop article s is the same as that of bulk NiO, except that the unit 
cell is slightly enlarged jHUMI The spherical geometry of NiO nanoparticle consists 
of circles stacked with decreasing radius on both sides of the equatorial great circle. 
These circles are circular cross-sections of (111) planes of NiO. The lattice sites 
in these circular planes are arranged in a triangular lattice structure. We show a 
part of a (111) plane in Fig. [2] The separation between two neighboring planes is 
8 = 2a/v / 3, where a is the triangular lattice parameter which is related to the cubic 
lattice parameter ao as a — ao/^/2. These circular planes in a NiO nanoparticle are 
stacked in a sequence A-B-C-A-B-C • • • as shown in Fig.[TJ where A, B, C planes are 
distinguished from each other by a shift of their centers from the origin. We label 
0,1,2 to the lattice points in the successive planes A, B, C. The position vectors 
of the center of planes A, B, C can be given (using cartesian unit vectors i and j) 
as r = 0, fx = a(| + j^i), anc ^ ^2 = a (§ — 57f^ res P ec ^i vei y- W e can write the 
three-dimensional position vector of the lattice sites in I th plane labelled by integers 
m, n and I (using the cartesian unit vectors k) as 

rw = {m + -)ai H —aj + l5k. (2) 
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If No, Ni and N% are the total number of lattice sites, counting from the planes of 
type A, B and C respectively, then the total number of lattice sites within a sphere 
of radius R will be given as, 

2 

N spheie (R) = ^N I (R), (3) 

7=0 

where, 

Ni(R)= ]T Q(R-\fmni~ri\), (4) 

r rn.nl 
l = I(Mod3) 

8 represents Heaviside step function — and 77, (J = 0,1,2), has been discussed 
above. We transform the above equation using the Poission sum formula ^^Sl as 

P-,q,w 

x Q(R-\? -fj|)dV. (5) 
Thus, the total number of spins can be written as 

16tt fR\ 3 ^ _ s ^ 31 (GR) 

0(1 



J EE-(w)s^ (6) 

1=° {6} 



where G = -42- 

V3a 



is a three dimensional reciprocal lattice 



vector, labeled by three integers p, q and w, and Ri = fi + /<5/e. ji is a spherical 
Bessel function of order one 

Due to the oscillatory behavior of the Bessel function, N sp here{R) varies non- 
monotonically with the particle size R, and the wavelength of oscillations goes as 
1/G. Thus, the longest wavelength mode G = in the above gives the smooth con- 
tribution as Adhere = T 71 " (ao ) ' wn ^ e the t erms with G^O represent oscillatory 

fluctuations. The fluctuation Nj p ^ere m ^ ne total number of spins in the sphere of 
radius R can be obtained from 



AT- ArBulk I iiifiuC jrj\ 

JVsphcrc — -<V sphcrc "f JV sphere . yt ) 

From Eq. (J6j> , and the asymptotic behavior of the spherical Bessel function, ji(x) ~ 
1/x, we can see that the amplitude of oscillatory fluctuations varies as, 

^sphere ~ R - ( 8 ) 

Hence the next to leading order term in the number of spins within a sphere goes 
as R rather than R 2 , which we could have expected from a random- walk argument, 
viz. the amplitude of the fluctuations is proportional to the square root of the the 
number of points on the boundary, here the spherical surface. The amplitude of 
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Fig. 4. The total magnetic moment M sp ^ ere (R) of the Neel-state for the sphere of radius R. The 
solid line is Af sp i lcrc (i?) obtained from analytical expression Eq. fl I L'li for the terms up to \G\ = 400 
whereas the dotted line shows the numerical counting using Eq. and Eq. JToL 

fluctuation varies as 1/G 2 . Thus larger G values give smaller contribution to the 
amplitude of fluctuation. Hence, we can get a good approximation to the total 
number of spins by retaining only few terms in the sum in Eq. (|6|). In Fig. [3j we 
have plotted the fluctuation in the number of spins as obtained from Eq. ([7|) for the 
terms up to \G\ = 400 along with that obtained from the exact numerical counting 
of spins. 

3.1. Magnetization fluctuations of the Neel state 

The bulk two-sublattice Neel magnetic structure of NiO requires us to assign all the 
spins in a circular plane of the NiO nanoparticle to be either +1 or -1. Thus we assign 
the circular planes to be +1 and -1 alternately, corresponding to ferromagnetic 
sheets of spins with alternating polarization stacked along (111) direction in the 
FCC lattice. Since, we have three different types of circular planes, and the circular 
planes are stacked as A-B-C-A-B-C • • • , as shown in Fig.[TJ all the spins in each of A 
(or B or C) type of planes will have either +1 or -1 value depending on the location 
of the plane along the stacking direction. The total magnetic moment M sp h e re(-R) 
of NiO spherical particles of size R can be found by summing the magnetic moment 
of all the circular planes. Following Hutchings et a/^2l, we assume that each Ni 2+ 
spin has a magnetic moment of 2/iB- Thus, we can write the total magnetic moment 
for the spherical particle as 

2 

-^sphere 
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where, 

M I (R)= (-lfQ(R-\rmnl-ri\). (10) 



l = I(Mod3) 



Applying the Poisson sum formula and proceeding analogously as we did above, 

M/(i?) = 37!w£ } cos(5 - r/)cos l^J 

x J e i8 p e^Q{R-r')<fr', (11) 

we evaluate Mj(R) and thus the total magnetic moment in the units of Bohr mag- 
neton (/zb) as 

. , 327T _R 2 ^ — ^ ^ — / _» x 

M sphcIC (R) = -^-—2^ cos (g- r i) 

a ° 1=0 {g,w} 

x cos ( ( G ? + 37r 2 + 2 ^ nG 



x jj ji? (g 2 + 3tt 2 + 2V3^G Z ) 5 1 , (12) 

where G = g + G z k and G z = §fw. The total magnetic moment Msp^ere displays 
oscillations as a function of the particle size, and the wavelength of oscillations goes 
as 1/ (G 2 + 37r 2 + 2V37rG z ) 2 . Unlike N sphcl . c , which had a smooth part (G = 0) 
and oscillatory terms (G =^ 0) (see Eq. ©), all the terms in the above Eq. (|12[) 
for the total magnetic moment display oscillations. In fact, all the terms have a 
similar asymptotic behavior. Using the asymptotic behavior of the Bessel function, 
the amplitude of the fluctuations in M sp i, ere can be shown to vary as, 

-^sphere ~ R- (13) 

The contribution from the longest wavelength mode G = to the magnetic moment 
can be written as 

(14 » 

The terms with (G ^ 0) represent the fluctuations in the total magnetic moment 
on various length scales. The magnetic moment, as obtained from Eq. (fT2|) for the 
terms up to |G| = 400 and as obtained from exact numerical counting is plotted 
with particle size in Fig. HI We find that the net magnetic moment is not as large 
as seen from experiments. For example, for the particles of diameter 3nm, we find 
the magnetic moment to be 26/ie, which is too small compared to experimental 
value 500/xb Also, for the particles of diameter 15nm, we find the magnetic 
moment to be 112/^b, whereas experimental investigation reports 700/xb Thus, 
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a net magnetic moment due to Neel-state does not quantify the large magnetic 
moment experimentally observed in NiO nanoparticles. 

In a realistic situation the shape of NiO nanoparticles is not a perfectly spher- 
ical. An oblate spheroid (with hi gh o blateness) or a platelet-shaped geometries 
has been reported by experiments! ^*^ Thus we again calculate magnetization for 
an spheroidal geometry of nanoparticles. The result of calculations for an oblate 
spheroid with polar radius R/e and equatorial radius R, shown in the Appendix, is 



MspheroidC-R) = -^-—3^2 X! COS (S-^) 
6 0,0 1=0 {g,w} 

{2Iirw\ 3tt 2 n y/3irG t 
x cos — — G + —tt- + 2- 



1 



2 




xji<R \G 2 + ^ r + - - }, (15) 



and the magnetic moment for the longest wavelength mode G = can be written 

as 





1 V^TT 














ao 1 



M 3=a =^=(-) . ! 1 { ' I (l'il 

Using this geometry, we find that the net magnetic moment value increases a little 
bit but again it is not comparable to experimental values, for 3nm particles, the 
magnetic moment improves to 34/ib for e = 2; 40/is for e = 3; and 70/iB for 
e = 4. These values are still very small as compared to experimental value 500/ie 
Similarly, for the particles of diameter 15nm, wc find the magnetic moment to 
be 153/zb for e = 2, 258^b f° r e = 3, and 308/ib for e = 4. Thus we see that 
increasing the oblateness increases the net magnetic moment of the nanoparticle. 
But the experimental value is still beyond our reach within the present model. In 
order to improve the model, we need to invoke a different ordering for surface spins 
than the bulk Neel-state ordering. We analyze the effects of roughness on the surface 
of nanoparticle and introduce a surface anisotropy and a ferromagnetic exchange 
interaction term for the surface spins in the Hamiltonian. The modified model is 
discussed in the following section and a variational approach is used to find the 
optimal thickness of surface roughness. 



4. Surface effects and variational approach 

The surface effects dominate the magnetic properties of nanoparticlesPEEEll 
breakdown of the dominant next-nearest neighbor antiferromagnetic interaction on 
the surface of the nanoparticle leads to uncompensated spins. These uncompensated 
spins play a vital role in determining the magnetic behavior of NiO nanoparticles. 
The magnetization reversal study for antiferromagnetic nanoparticles by Zianni et 
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al using Monte Carlo simulation — also reveals a distinct magnetic role of surface 
and core spins. The broken bonds and defects at the surface layer gives rise to high 
magnetic response of the disordered surface spins than the core spins j22| 23 | 24 | 25 | 
This surface effect can be incorporated in two different ways; viz. an ordered surface 
with surface spins ordering owing to internal field due to core spins, and a disordered 
surface with spins oriented randomly at the surface which cancel out thus giving 
zero contribution to the total magnetic moment. Due to the surface roughness, the 
uncompensated surface spins can be more easily polarized by a small magnetic field. 

For an ordered surface, we take that the spins inside a surface roughness shell of 
thickness A are aligned by a field due to a net core magnetic moment though small. 
This would enhance the net magnetic moment of nanoparticles. In this scenario, the 
core spins within a sphere of size R — A have the bulk antifcrromagnetic structure, 
carrying a magnetic moment of order R, as we calculated in Sect. [3j and the spins 
within a shell of size A are all polarized, carrying a magnetic moment of order R 2 . 

For a nanoparticle of radius R, we can write the total magnetic moment 

cicrcd 

(R, A) as 

M ordere d (R, A) = \M sphcrc (R - A)| + 2/i B (^sphere (R) 

- N spheIe (R- A)). (17) 

In the above, the first term is due to the core Neel-state magnetic moment, and 
the second term represents the surface roughness effect. Here, the core spins within 
a sphere of radius R — A have the bulk magnetic structure. The spins within the 
shell of thickness A are aligned, each spin contributing a magnetic moment of 2//b- 
Since M sp h C rc in the Eq. (JTTj) goes as R while surface roughness terms as a whole is 
proportional to i? 2 , the total magnetic moment has a leading term going as R 2 , if 
the shell thickness is independent of the size. 

We take a variational approach to establish the behavior of A with the particle 
size. In the variational approach, we modify the Hamiltonian (Eq. (fTJ)) for nanopar- 
ticle, by including a ferromagnetic exchange interaction term —Js^2uj\ Si-Sj with 
coefficient J s for spins lying in the roughness shell A and a ferromagnetic exchange 
interaction term like —J cs ^2 iecorc Si.Sj with coefficient J cs between core and sur- 

j'G surface 

face spins lying at the interface of core and surface. Moreover we introduce surface 
anisotropy term — K s X^(A-^) 2 with coefficient K s for the surface spins which pre- 
vails over the uniaxial core anisotropy with coefficient K\ . Hence we can ignore the 
core anisotropy term. Taking all the interactions into account exactly, we can write 
energy of the system as 

73(A) = {fi- a)N sphcl . c (R - A) - pN sphclc (R) 

- E ^>< ( 18 ) 

z£corc 
j£ surface 

where a — z c + 2J + 2 j"" and j3 — z s -j^ — h -j 3 . z c and z s are the coordination 
numbers of core spins and surface spins. We use the value of J cs same as that of 
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Fig. 5. The surface roughness thickness A is plotted with the particle size. Dotted line shows the 
optimized A obtained from Eq. I|18|l . The fluctuations in A/R show an envelope decay, A/R ~ 
.26R~ 6 / 5 , shown as a solid line. 



J nnn ■ In the above equation, we find that the hrst term and the last term explicitly 
contain the variational parameter A. The energy is minimized with respect to A to 
get the optimal thickness of the surface roughness shell. We choose the parameters 
to be a = 6.24ao/i? and f3 = 6ao/R, such that we have a nontrivial A and the total 
magnetic moment per particle for an assembly of particles with lognormal distribu- 
tion is close to the experimental value for the corresponding sizes. This indicates 
a modification in the core interaction strength in nanoparticles, which signifies a 
deviation from a two-sublattice ordering. For large enough sizes, A becomes zero 
and only contribution to the net magnetic moment is the first term in Eq. (|17l) . 

In Fig.[5l the optimal A/R which minimizes energy E in Eq. (| 18[) is plotted with 
particle size R. We find that A/R shows an oscillatory behavior and the amplitude 
of oscillations is decreasing with increasing particle size. The best fit of the curve 
shows that A ~ .26.R -1 / 5 . In Fig. El we have plotted net magnetic moment for 
corresponding particle sizes using Eq. (|17j) . The total magnetic moment displays 
size-dependent fluctuations, whose amplitude shows a peak at R ~ 12a . Further 
increasing the size of the particle lowers the magnetization which occurs due to the 
lowering of the surface roughness effect i.e., less availability of ordered spins near 
surface. A similar behavior has been observed in experiments by Yi et al ^ where 
the magnetization of NiO powder increases with the annealing temperature (grain 
size) and shows a peak at an annealing temperature 170°C. Annealing at higher 
temperature than 170°C leads to a lower magnetization. 

In the case of a disordered surface, we take a disordered spin structure inside 
the shell of thickness A. In this situation, the net contribution from the disordered 
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Fig. 6. The net magnetic moment as a function of particle size, including the surface roughness 
and anisotropy, calculated from Eq. l|17|l . For each size an optimized surface roughness thickness 
A from Eq. H18H is used. 

surface will be zero and magnetic moment will solely be due to the core spins inside 
the sphere of radius R — A. 

^disordered = Msphere(-R — A). (19) 

Since M sp hore shows nonmonotonic oscillatory dependence on particle size R, 
-^disordered should also show the same behavior with particle size. In Fig.[7][a) we plot 
coarse-grained -Mdi SO rdered with particle size R using spherical as well as spheroidal 
particles. Averaging is done over a window size of 2.8cto. In the same curve we have 
also shown ordered surface case. The net magnetic moment increases with increas- 
ing oblateness of the spheroidal nanoparticles. In the same Fig. [7f a) we have shown 
ordered surface case for both geometries, spherical(e = 1) and spheroidal with e = 4. 
In both the geometries, the coarse-grained magnetic moment shows nonmonotonic 
oscillating behavior as a function of the size. Though the peak value in spheroid 
case is more than the sphere, we see a sharp rise and fall of magnetic moment just 
before and just after the peak (w 9ao) in spheroid case which is different than sphere 
where the magnetic moment slowly increases with particle size, a plateau is seen for 
intermediate range, and for bigger sizes it starts decreasing with size. 

The systems of magnetic nanoparticles in experimental studies are in general 
polydisperse. The shape and size of the particl es a re not well known but the particle 
size distribution is often found to be lognormal^l We consider the system consisting 
of lognormally distributed, widely dispersed nanoparticles, hence non interacting 
among each other. The weight of a given size (radius R) of nanoparticles is given 
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by a lognormal distribution, 



P(R/a ) 



1 



e 



(l n (R/a )-^) 2 



(20) 



a{R/a )V^ 



The characteristic parameters of the distribution are chosen to be /i = ln(R/ao) and 
(j = 0.5. Using this distribution, the magnetic moment for ordered-surface spherical 
nanoparticle of diameter 3nm works out to be 510/ib which is quite close to the 
experimental value 500/xb For the ordered-surface oblate spheroidal nanoparticle 
(oblateness e = 4) of same size this value is 297/ib, which is a little less than the 
experimental value. Using the same size distribution for disordered-surface spherical 
nanoparticle and spheroidal nanoparticle, we find net magnetic moment 18/iB and 
40/ib respectively, which is very small as compared to experimental value. Similarly, 
for a distribution of particles with mean size 15nm, our calculation of net mag- 
netic moment for ordered-surface spherical particle is 820//B, and for order-surface 
spheroidal particle is 703/iB- For disordered case, calculated value for spherical and 
spheroidal geometries are 68/ib and 122/^b, respectively. Thus for 15nm size parti- 
cles, magnetic moment calculations using ordered surface and spheroidal geometry 
is very close to the experimental value 700/LtB-^The experimental value of magnetic 
moment for 8.5nm particle size^^Sl a t \QK (700/ib) compares quite well with our 
calculated value 802/ib with ordered-surface and spherical geometry. However, mag- 
netic moment value using particle with ordered-surface and spheroidal geometry is 
550/ib, a little less than the experimental value. In our model, the total magnetic 
moment per particle may slightly increase or decrease depending on a distribution 
of particle sizes and oblateness in the shape of particles. 

In Fig. [2Kb) we have shown the total magnetic moment averaged over a log- 
normal distribution defined above as a function of the size of the particle. For an 
ordered case starting from very small sizes, we find that the magnetic moment 
increases with increasing the particle size and attains a peak, and further increas- 
ing the particle size decreases the net magnetic moment. The averaged magnetic 
moment for a spheroid particle with e — 4 has also been shown in the Fig. [7jb). 
For a spherical particle magnetic moment has peak value around 12ao whereas for 
spheroid case it is around 9ao-As we can see from Fig.[5l the surface roughness effect 
is stronger only for particles of intermediate sizes, which confirms the greater role 
of surface roughness for these sizes. But for the nanoparticles of sizes greater than 
12eto, the surface roughness shell A becomes very small and a net magnetic moment 
arises largely due to the uncompensation of bulk N eel-state ordering which will tend 
to zero for large enough sizes showing the bulk character. Thus Fig. E^b) reflects a 
trend for ordered surface particles, where the net magnetic moment is very small for 
smallest size particles. Increasing the size, magnetic moment increases and reaches 
a maximum due to the greater role of surface, and again decreases towards the bulk 
value. In the same figure FigJT^b), we have shown total magnetic moment for a 
disordered case averaged over the same distribution as discussed above. The mag- 
netic moment in disordered case is showing an increasing trend with particle size, 
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Fig. 7. (a) Coarse-grained magnetic moment versus particle size where coarse-graining is done 
over a window of size 2.8arj. (b) Averaged magnetic moment versus size where averaging is done 
over a lognormal size distribution, shown in Eq. (1201) . with a width of the distribution a = 0.5. 
In both figures (a) and (b), two different cases ordered and disordered surfaces has been shown. 
For both cases two different geometries of nanopartcles, a sphere(e = 1) and a spheroid with e = 4 
has been shown. An increase in oblateness e results an enhanced magnetic moment in disordered 
case, but for ordered surface, magnetic moment for spheroidal particle is less than that of spherical 
particle. 

but the value is small as compared to ordered case. The magnetic moment values 
calculated here depends on the model parameters the core interaction parameter a, 
and the surface roughness parameter ft, and the width of the size distribution a. 
Adjusting these parameters for a given size distribution and varying oblateness leads 
to the net magnetic moment value comparable to the experimental values. Tuning 
of the parameters a and f3 are directly related to the multi-sublattice ordering for 
nanoparticles as predicted by Kodama et al\^ 

5. Conclusions 

We have investigated finite-size and surface roughness effects in NiO nanoparti- 
cles. We have found that the net magnetic moment due to finite-size fluctuations is 
nonmonotonic, oscillatory and proportional to the particle size R, hence magneti- 
zation goes as 1/R 2 . The geometry of the particle also plays an important role in 
net magnetic moment. An oblate spheroid shape particle shows an increase in net 
magnetic moment by increasing oblateness of the particle. The experimental mag- 
netic moments for various sizes are quite large compared to the magnetic moments 
that arise as a finite-size fluctuation. The surface effects become very important in 
nanoparticles. We have incorporated surface effects in two different ways; an ordered 
surface where all the spins lying in the surface roughness shell are aligned due to 
internal field from the core spins, and a disordered surface where spins are randomly 
oriented in the surface roughness shell. Due to roughness of the surface and struc- 
tural disorders, the uncompensated surface spins can be more easily deviated from 
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the antiferromagnetic alignment by a magnetic field. We have introduced a surface 
anisotropy term and a ferromagnetic exchange interaction term for the surface spins, 
and a ferromagnetic exchange interaction term between core and surface spins lying 
at the interface of core and surface along with the bulk model in the Hamiltonian. 
A variational approach has been taken to find the dependence of the shell thickness 
on the size of particle. We have found that for nontrivial values of A, the core in- 
teraction strength is modified which shows a signature of multi-sublattice ordering 
rather than two-sublattice ordering for smaller sizes. A is showing size dependent 
fluctuations, with an envelope decay A ~ i? -1 / 5 . We have shown that the total 
magnetic moment calculated with ordered as well as disordered surfaces displays 
size dependent fluctuations. For an ordered surface case, smoothening these fluctu- 
ations by a window-averaging, using a lognormal size distribution of nanoparticles, 
results a magnetic moment per particle which is very close to observed experimen- 
tal values of various sizes. We have also found that due to surface roughness effect, 
the net magnetic moment shows a trend where magnetic moment is very small for 
smallest size particles. Increasing the size, magnetic moment increases and reaches 
a maximum at R ~ 12ao (depending upon the distribution of sizes and oblateness 
of the particles), and again decreases towards the bulk value. 
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Appendix A. Spheroidal nanoparticles 

In this appendix, we calculate the total number of spins and net magnetic moment 
of a nanoparticle using spheroid geometry with polar radius R/e and equatorial 
area R. We can write the three-dimensional position vector of the lattice sites in I th 
plane labelled by integers to, n and / is (using the cartesian unit vectors k) as 

f m ni = (m + -)ai H —aj + lSk. (A.l) 

The spheroidal geometry implies that only those spins whose locations satisfy xL , + 

V 2 mnl + ^^nnl < ^ \ sh ° uld be Counted. 

If No, Ni and N2 are the total number of lattice sites, counting from the planes 
of type A, B and C respectively following section [31 then the total number of lattice 
sites within a spheroid of equatorial radius R will be given as, 

2 

-^spheroid (R)=Y,Ni(R), (A.2) 
7=0 
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where, 

N I (R) = Yl 6(_R 2 -( SrallI -.i;/) 2 -(fc„i-w) 2 -e 2 i 1 ), (A.3) 

mnl 

l=I(Mod3) 

and xo = 0, yo = 0; x\ — a/2, y\ — a/2V3; and X2 — a/2 , y? = — a/2v3. We 
transform the above equation using the Poission sum formula l^SI as 



Ti(x'p+y'q+z'w) 



X © f A - \l {x'mni - %i) 2 - (y' mn i ~ Vi) 2 ~ e ^'Li) dx' dy' dz' . (A.4) 
Thus, the total number of spins can be written as 

N sphermd (R) = ^ ^g^, (A.5) 

6 \ a °^ 1=0 0} 



where G = -4§- 



is a three dimensional reciprocal lattice 



vector, labeled by three integers p, q and w, and Ri — rj + IeSk. j\ is a spherical 
Bessel function of order one ^1 

We can write the total magnetic moment for the ellipsoidal particle as 

2 

-^spheroid (R)=2hb^2Mi{R), (A.6) 

7=0 

where, 

M I (R)= J2 (-lM(^), (A.7) 

r mn ; 
i = /(Mod3) 



and Ri = \/ R 2 — F<5 2 e 2 . Applying the Poisson sum formula and proceeding analo- 
gously as we did above, we calculate Mj(R) and thus the total magnetic moment 
in the units of Bohr magneton (/xb) as 

M sph eroid(-R) = -^-—s'yi E COs(5.f» 

3e ao' 3 f-^ ^ 

1=0 {g,iu} 



X COS 



xj 1 {R[& + ^ + --^l) >, (A. 



, 3tt 2 , 2\/37rG 

n, i l 

where G = g + G z k and G 2 = ^rjw- 
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